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Abstract Minor bodies of the solar system can be used to measure the spectrum 
of the Sun as a star by observing sunlight reflected by their surfaces. To perform 
an accurate measurement of the radial velocity of the Sun as a star by this method, 
it is necessary to take into account the Doppler shifts introduced by the motion of 
the reflecting body. Here we discuss the effect of its rotation. It gives a vanishing 
contribution only when the inclinations of the body rotation axis to the directions 
of the Sun and of the Earth observer are the same. When this is not the case, the 
perturbation of the radial velocity does not vanish and can reach up to ~ 2.4 m/s 
for an asteroid such as 2 Pallas that has an inclination of the spin axis to the 
plane of the ecliptic of ~ 30°. We introduce a geometric model to compute the 
perturbation in the case of a uniformly reflecting body of spherical or triaxial 
ellipsoidal shape and provide general results to easily estimate the magnitude of 
the effect. 

Keywords Techniques: radial velocities - methods: data analysis - Sun: general 
- Sun: photosphere - minor planets, asteroids: general. 


1 Introduction 


Obtaining a spectrum of the Sun integrated over its disc, i.e., directly comparable 
with stellar observations, is not an easy task. Generally, the spectrum reflected by 
a minor body of the solar system or by one of the Galileian satellites has been used 
as a proxy for the spectrum of the Sun as a star (cf. |Molaro Gentnrionl 2011). 
An accurate measurement of the wavelength of a spectral line in a reference frame 
at rest with respect to the barycentre of the Sun, or an accurate measurement of 
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the radial velocity (hereafter RV) of the Sun as a star, require a correction for the 
Doppler shift produced by the motion of the reflecting body (cf. [Molaro Monai] 
2Q12[ Molaro et al.[ |2Q13| ). The effect of the orbital motion with respect to the 


barycentre of the Sun and to the observer on the Earth can be corrected using the 
NASA Horizon ephemeride^ but the effect of the axial rotation of the body must 
also be taken into account when a precision of the order of 0.1 — 1 m/s is required. 
This is the case of the observations of the Sun as a star performed to understand 
the impact of solar convection and magnetic activity on its disc-integrated RV. 
These investigations are of fundamental importance to understand similar effects 


in distant solar-like stars that are searched for Earth-like planets (e.g., Lanza et 
2Q11[ Dumusque et al.[ 2012). 


al 


This work is dedicated to a precise computation of the effect of the axial 
rotation of a reflecting body - hereinafter indicated, for simplicity, as an asteroid 
- on the solar radial velocity. Eor simplicity, we shall assume that the body has a 
spherical surface in Sect. [2^ while the case of a triaxial ellipsoidal shape will be 
considered in Sect. 12.41 

We consider a body of uniform albedo. Patches with a different albedo on the 
surface of a rotating asteroid may produce distortions of the line profiles in the 
reflected spectrum that are modulated with the rotation period of the body. The 
origin of the line distortions is the different level of the continuum in the spectrum 
reflected from an albedo inhomogeneity. Eor example, in the case of a dark spot, 
the local lower continuum level produces a bump in the spectrum integrated over 
the disc of the asteroid as in the case of a spotted star (cf. Eig. 1 in Vogt & 


Penrod 1983). This affects the measured solar RV with a periodic perturbation 


having the same period of the rotation of the asteroicQ It is possible to correct for 
this effect by fitting a sinusoid and its harmonics to the RV time series with their 
fundamental period equal to the rotation period of the asteroid. The coefficients 
of this Eourier series will be a slowly varying function of the viewing angle of the 
asteroid spin axis that implies that this method can be applied only for time series 


that are much shorter than the asteroid and the Earth orbital periods (Haywood 


et al.| 2015). On the other hand, the effect that we investigate in the present work 
is not modulated with the rotation of the asteroid, but varies slowly with the 
relative position of the body and of the Earth along their orbits. Therefore, it is 
not possible to correct for this effect by the simple method applicable in the case 
of the modulation arising from the albedo inhomogeneities. 


2 Model 

2.1 Reference frame 

We consider a reference frame with the origin at the barycentre T of the asteroid. 
Its spin axis points in the direction {is,j3s), where £s is the ecliptic longitude and 

^ http://ssd.jpl.nasa.gov/?horizons 

^ In principle, this effect can be used to measure the rotation period of an asteroid if the 
amplitude of the modulation is comparable or larger than the intrinsic radial velocity variations 
of the Sun on typical rotational timescales, i.e., a few hours or days. Those i ntrinsic variations 
are dominated by surface convection and have amplitudes of a few m/s (cf. [Dumusque et al.j 
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13s the ecliptic latitude of the projection of its North pole onto the celestial sphere 
as seen from T. The ecliptic longitude of the Sun and of the (geocentric) observer, 
as seen from T at the epoch of sunlight reflectiorj^ are indicated with Iq and 
respectively. The inclinations of the asteroid spin axis to the direction of the Sun 
and to the observer lying on the plane of the ecliptic, i© and Z 0 , respectively, are 
given by: 


cosz© = cos (-£0 —-£s) cos/3s, ( 1 ) 

cos*© = cos(£© -£s) COSTS' ( 2 ) 

The ecliptic longitude of the Sun and of the geocentric observer as seen from T 
can be obtained from the longitude of the asteroid as seen from the barycentre 
of the Sun -£o© and of the Earth ^o©, respectively, given by the NASA Horizon 
ephemerides, as: Iq = tt — -^o© and = tt — -£o©- From Eqs. Q and § , we see 
that both the angles z© and z© vary between /3s and tv — /3s as the asteroid revolves 
around the Sun. 


2.2 Radial velocity variation induced by the asteroid rotation 


We consider a Cartesian reference frame with the origin at T, the z-axis along 
the spin axis of the asteroid and the x and y axes in the equatorial plane of 
the asteroid. In addition to the Cartesian coordinates, we consider also spherical 
coordinates (cf. Fig. [^. The origin of the longitude is chosen in such a way that 
the longitude of the centre of the Sun is — o/ 2 , while that of the observer is a/2 at 
the time of sunlight reflection (cf. Fig.[^ - a similar approach has been introduced 
to model asteroid photometric variations (cf., [Harris et al.[ 1984). Therefore, the 
unit vectors from T to the centre of the Sun S and to the geocentric observer O 
have Cartesian components: 


TO = (sinz© cos(a/2), sinz© sin(o/2), cosz©), 
TS = (sinz© cos(a/2), — sinz© sin(o/2), cosz©). 


( 3 ) 

( 4 ) 


A generic point on the surface of the asteroid has Cartesian coordinates P = 
R(sin^cos A,sin^sin A, cos^), where 0 < ^ < tt is its colatitude measured from the 
North pole, 0 < A < 27r its longitude, and R the radius of the body assumed to be 
spherical (see Fig.[^. Let us consider the angles '0© between the normal in P and 
the direction to the observer TO, and V^© between the normal and the direction 
to the Sun TS. By performing the scalar products between the unit vectors OP 
and TO, and OP and TS, the Cartesian components of which are given above, we 
find: 


/z© = COS'?/;© = sinz© sin^cos(A — a/2) + cos^cos^, (5) 

/z© = COS'?/;© = sinz© sin^cos(A P a/2) + cos^cosz©. (6) 

^ For simplicity, we can neglect the delays due to the time taken by sunlight to travel from 
the Sun to the asteroid and be reflected to the observer because the intervening variations in 
the positions of the bodies are too small to signiflcantly affect the angles entering in our model 
(see below). 
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The rotation velocity of the point P is found by differentiating its position vector 
as a function of the time because the longitude of P increases steady in our fixed 
reference frame owing to the rotation of the asteroid. Introducing the equatorial 
rotation velocity f/gq = (27r/Prot)^, where Prot is the rotation period, we find: 

V(P) = i4q(— sin sin A, sin ^ cos A, 0). (7) 

The measured RV is a weighted average over the illuminated portion of the as¬ 
teroid’s disc, where the weight of each disc element is proportional to the flux 
received from it. Therefore, we need to define the limits of the visible disc from 
which a non-zero flux is received by the observer. The limb of the visible disc is 
defined by the condition = 0 that gives the longitude limits A 0 (^)(^) of the disc 
for a given colatitude 0 as: 


cos(A 0 ('^) — o/ 2 ) = — cot % cot 0, 2 = 1 , 2 . ( 8 ) 

On the other hand, the limb of the illuminated hemisphere of the asteroid is 
defined by the condition: fiQ = 0. Indicating with XQ(^k) f^ie longitude limits of the 
illuminated hemisphere at a given colatitude 6 >, we find: 

cos(A 0 (/j.) + o/2) = — cot 20 cot 0, k = 1, 2. (9) 

We introduce the angles 70 = arccos(— cot 20 cos 0) and 7 © = arccos(— cot 2 © cot 0). 
For a given colatitude 0 on the disc of the asteroid, the visible and illuminated lon¬ 
gitude range [Ai(^), A 2 (^)] is given by (cf. Fig. where we show the case of the 
equatorial plane, i.e., = 7r/2): 

[Ai, A2] = [max(-70 + a/ 2 , -7© - a/ 2 ),min (70 + 0 / 2 , 7 © - a/ 2 )]. (10) 


A given point P on the illuminated and visible disc of the asteroid receives the 
light coming from the Sun with a Doppler shift corresponding to its radial velocity 
with respect to the barycentre of the Sun, FR(rot)i(^)? reflects the spectrum 
towards the observer without any further shift in its rest frame. For simplicity, we 
consider here only the component of the radial velocity produced by the rotation 
of the body, that we indicate with the subscript (rot). The observer on the Earth 
is moving with respect to the reflecting point with a radial velocity FR(rot)( 2)5 
also due to the rotation of the body, that introduces a further Doppler shift in 
the spectrum. In conclusion, the observed spectrum coming from P is Doppler 
shifted with a radial velocity corresponding to the sum of the two contributions, 

i-e-, y R(rot)(l){P) + ^(rot)( 2 )(^) (cf- 
2013 ). 

The perturbation of the measured RV produced by the rotation of the reflecting 
body is found by averaging the above Doppler shift over the surface of its visible 
and illuminated disc, i.e.: 


N ieuwenhuij zen 


1969 


Gjurchinovski 


2005 


ARV = 


lD,[yR(rot)(l)(P) + V^rot)(2){P)]dF(P) 

lo^dFiP) 


( 11 ) 


where dF{P) is the flux coming from the area element of the disc around the point 
P; the integration is extended over the illuminated portion of the disc Dj. 
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The flux coming from a given surface element at frequency v is given by: 

dF{P) = l{v,P)p.^{P)dA{P), (12) 

where l{y^ P) is the specific intensity of the element at u and dA{P) = sin 0 dOdX 
the area of the surface element. 

Since the asteroid is illuminated by the Sun, the intensity reflected at a given 
point of its surface is given by Lambert’s law (e.g., |Kopal[ [1959 ): 

I{iy,P) = Io{iy)a{iy) cos'ipQ = Io{u)a{iy)fiQ, (13) 

where Io(iy) is the intensity for normal reflection and a{u) the albedo at frequency 
u. The total radial velocity difference produced at the point P by the rotation of 
the asteroid is the sum of the projections ldi(rot)(i) = ^R(rot)(2) = 

V(P) • TS] therefore, we obtain: 

^R(rot)(i) + '^R(rot)(2) = “Veq sin 6> [sin % sin(A - a/2) + sinz© sin(A + a/2)] . (14) 





Fig. 1 The adopted Cartesian reference frame with its origin at the barycentre T of the 
asteroid, the z-axis along its spin axis and the x-y plane in its equatorial plane. The spherical 
coordinates of a generic point P on its surface are also indicated - 9 colatitude from the North 
pole; A longitude; P' is the projection of P on the equatorial plane. 
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The integral in the numerator of Eq. (11) at a given frequency v becomes: 


/ [ 

JDy 


[^R(rot)(l)(^) + ^(rot)(2)] dF{P) — 


nTT PA 

= -Voc,h{v)a{v)R^ 

•Jo •J Ai 

+ sinz© sin(A + o/2)]/re//© , 


A2(0) 


(0) 


o^Asin^^ [sinz© sin(A — o/2) 


(15) 


where the limits of integration Ai and A2 have been specified above. In general, 
this integral can be evaluated only numerically owing to the non-closed expressions 
giving the integration limits Ai,2 as a function of 0. The total fiux F at frequency u 
received from the disc of the asteroid that appears in the denominator of Eq. (11), 
becomes: 

rA2(0) 

dA/Z0/z© sin^, (16) 


F = Io{v)a{v)R^ 
that can also be integrated numerically. 


r 

•Jo J Ai 


(0) 



Fig. 2 The sub-observer point O and the sub-solar point S in the adopted reference frame with 
the origin at the barycentre T of the asteroid. The initial meridian from which the longitude is 
measured is chosen so that the longitude of the points S and O be — 0/2 and o/2, respectively, 
where a is the angle between the directions to the observer and to the Sun with its vertex 
at the barycentre of the asteroid T (origin of the Cartesian reference frame). The inclinations 
of the TO and TS directions to the spin axis of the asteroid are indicated as z© and z©, 
respectively. 
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In conclusion, the radial velocity perturbation is given by: 


ARV = 

frT ^ [sin 20 sin(A — a/2) + sinzQ sin(A + a/2)]fi^fiQ dO dX 

= -( 17 ) 

fo I\ /e) Me M© sin 6» d6» dA, 

that is independent of the albedo and the specific intensity at the given frequency. 
Therefore, this RV shift can be applied to the whole spectrum. 


2.3 A particular case 

When the spin axis of the asteroid is orthogonal to the plane of the ecliptic, i.e., 
/3s = 7r/2, Eqs. 0 and 0 gives: 20 = 2© = 7r/2. A more general case arises when 
the longitudes of the Sun and of the observer are such that 20 = 2© = 2, where the 
common inclination 2 is arbitrary. Now, we shall consider this particular case. 



Fig. 3 Projections of the TO and TS directions onto the equatorial plane of the asteroid 
(6 = 7r/2). The extreme longitudes marking the limb of the observable disc (A0('x 2)) ^^iid 
of the illuminated hemisphere (Aq^^^ 2)) indicated, respectively. The limits of integration 
to be considered in the RV and flux calculations are Ai = max(A0i,A0i) = A01 and A2 = 
min(A02,A02) = A02, in this speciflc case. For 0 7^ 7r/2, the determination of the limits of 
integration is similar (see the text). 
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Comparing Eqs. (§ and we see that: 


^©(i) = -^0(/c) with i ^ k. 


(18) 


For a given colatitude 6^ the longitude of the illuminated portion of the visible disc 
ranges from Ai(^) = to X 2 ( 0 ) = X( 0 )q(^ 2 ) when the limb of the visible disc 

at A(^)0('i) is illuminated (cf. Fig.j^, or from Ai(^) = A(^)q(i) to X 2 { 0 ) = X{ 0 )^(^ 2 ) 
when the limb at A(^)0(2) is illuminated. 

In the particular case that we are considering, the integral at the numerator 
of Eq. becomes: 


/1 

JDi 


FR(rot)(l)(-P) + ©ot)(2)] dF{P) = 


= —VeqIo{^')a{iy)R^smi / sin^ 0 [sin(A — a/2) + sin(A + a/2)]/X0/xo d0dA (19) 


JDi 


P7T rX2{0) 

=—VeqIo{iy)a{u)R^ sini / dO dA sin^ ^ [sin(A — a/2) + sin(A + a/2)]/X0/i0 , 

Jo Jxi{e) 

where Ai and A2 were specified above and are equal to zero when the colatitude 0 
corresponds to a point outside the visible disc of the asteroid. 


The evaluation of the integral (19) can be made by considering the symmetry 
of the integrand function with respect to a change in the sign of A. Considering 
Eqs. and the product is invariant under a change in the sign of A. 

On the other hand, the same transformation changes the sign of [sin(A — a/2) + 
sin(A + a/2)] thus making the integrand in Eq. (19) antisymmetric with respect 
to the transformation A ^ — A. When we perform this change of variable in the 
integral, the limits of integration changes according to Eq. (18), say, becomes 

A0(2) and A0(2) becomes A0(i). Together with the sign change in the differential 
of the integration variable dA, this results in no change in the limits of integration. 


This implies that the integral (19) vanishes. In other words, by integrating the 


Doppler shift over the visible disc of the asteroid, we find that the net change in 
the solar RV is zero when = ^o- 


2.4 An asteroid of ellipsoidal shape 


The model in Sect. |2.2| can be extended to the case of an ellipsoidal asteroid whose 
equation referred to the Cartesian reference frame with origin at T is: 


f{x,y,z) 



( 20 ) 


where x,y,z are the Cartesian coordinates of a point P on the surface of the body 
and tte, be, and Ce are its semiaxes (cf. Fig. [^. In the adopted spherical polar 
reference frame, the Cartesian coordinates of P can be expressed as a function of 
its spherical coordinates as: P = (aesindcosA,5esindsinA,CeCosd). Using these 
parametric equations for the ellipsoidal surface, we can compute the Cartesian 
components of the elementary vectors tg and tx that are tangent to the surface 
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at the point P in the colatitude and in the longitude directions, respectively. By 
differentiating the equations, we find: 

t 0 = (tte COS A COS 0, be siu A cos 6>, — Ce sin 0)d0 and (21) 

ty = (—aesinAsin6>,5eCos Asin6>,0)dA. (22) 


In general, ty and tg are not perpendicular to each other, except when a© = be. 
The elementary surface area at the point P can be obtained as the modulus of 
the cross product of the two tangent vectors, i.e., dA = jt^ x ty| that can be easily 
computed from their Cartesian components. 

The unit normal h at P is parallel to Vf{x,y, z). Expressing the Cartesian 
components of the gradient by means of the spherical coordinates: 


HP) = i ( 


sin 0 cos A sin 0 sin A cos 0 


ae 


Ce 




(23) 



Fig. 4 Illustration of the ellipsoidal model adopted to compute radial velocity variations in 
Sect. |2.4| The semiaxes of the ellipsoid ae, be, and Ce are labelled and are assumed to be 
parallel to the axes of the Cartesian reference frame. The spherical polar coordinates A and 9 
of a generic point P on the surface of the ellipsoid are indicated together with the local unit 
normal h and the unit vectors ty and tg tangent in P in the directions of increasing longitude 
and colatitude, respectively. P' is the projection of P onto the equatorial plane of the asteroid 
with its rotation axis assumed to coincide with the z-axis. 
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where 


G = 


A 


+ 




■ 2a I COs2 6» 

Sin 0 H-^—. 


(24) 


The expressions of the unit vectors TO and TS do not change for an ellipsoidal 
asteroid (see Eqs. [^and[^, while the projection factors = h-TO and jiQ = ri-TS 
can be obtained from the above components of the unit normal h at each surface 
point. 

The components of the rotation velocity of a given point P(A, 0) on the surface 
of the ellipsoid can be obtained by differentiating its longitude with respect to the 
time, as in the case of a spherical body. We obtain: 


V(P) 


—— (—cie sin A sin 6>, he cos A sin 0, 0), 

Prot 


(25) 


where Prot is its rotation period. This equation can be used to compute the com¬ 
ponents of the rotation velocity that appears in Eq. (11) by performing the scalar 
products as in Sect. |2.2[ 


The integrations in Eq. (11) can be performed numerically by dividing the 
ellipsoid into many surface elements and computing the contribution of each ele¬ 
ment to the radial velocity perturbation, i.e., to its numerator, and to the total 
flux, i.e, to its denominator. The expression for the elementary flux is not changed 
and is given by Eqs. (12) and (13). Einally, the integrals are obtained by summing 
the contributions of all the surface elements that are both illuminated by the Sun 
(M0 > 0) and visible from the observer (//0 > 0). 


3 Results 


Eor some bright asteroids, we list the mean radius, perihelion distance from the 
Sun, rotation period, equatorial rotation velocity, ecliptic longitude and latitude 
of the spin pole with respect to the J2000 reference frame, and the corresponding 
reference in Table Some of the rotation periods are taken from the NASA 
Horizon Ephemerides. Eor 20 Massalia there is an ambiguity in the longitude 
of the pole As, so both values are listed. In our case, the maximum value of 
Eeq = 92.27 m/s is achieved in the case of 1 Ceres. 

The angle a between the direction of the observer and that of the Sun at 
the time of light reflection is virtually identical to the Sun-Target-Observer angle 
(STO) as given by the Horizon ephemerides. The maximum value of a is estimated 
as Omax ~ arcsin(1.016/^), where q is the minimum heliocentric distance of the 
asteroid in AU. This conflguration corresponds to the Earth and the Sun seen in 
quadrature from the asteroid, while it is at the perihelion and the Earth at the 
aphelion. Using the data listed in Table we And Omax ~ 31°. The maximum 
values of ARV are listed in Table in the case of a spherical body. They were 
obtained by looking for the maximum radial velocity variation as given by Eq. (17) 
when the difference Iq — is is varied from 0° to 360° and the difference ^0 — = 

(i(7) — ^s)=t<^max because these two extreme values lead to the largest radial velocity 
perturbation. 

The ratio of the maximum RV perturbation to the equatorial velocity vs. the 
latitude of the spin pole is plotted in Eig. for different values of the angle o. 
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Asteroid 

q 

(AU) 

R 

(km) 

■Prot 

(hr) 

Veq 

(m/s) 

(deg) 

(deg) 

Reference 

1 Ceres 

2.55 

480 

9.075 

92.27 

346 

82 

1 

2 Pallas 

2.12 

275 

7.811 

63.65 

30 

-16 

2 

3 Juno 

1.98 

117 

7.210 

28.31 

108 

38 

3 

4 Vesta 

2.14 

265 

5.342 

86.53 

336 

63 

4 

7 Iris 

1.84 

100 

7.139 

24.44 

15 

25 

5 

20 Massalia 

2.06 

73 

8.098 

15.72 

31/208 

69 

3 


Table 1 Asteroid parameters. References: 1: 

Drummond et al. 

(20141; 2: 

Carry et al. 

(2010 

3:|Dotto et al.|(|1995|); 4:|Thomas et al.|(|1997 

1 ; 5:|de Pater et al.| (|1994|). 


considering a spherical body. This is an adequate assumption for large, bright as¬ 
teroids that deviate from a spherical shape by a negligible amount for our purpose. 
For instance, in the case of 1 Ceres, the surface has the shape of an oblate spheroid 
with semiaxes ae = be = 967 ± 10 km and Ce = 892 ± 10 km, i.e. an oblateness less 


than 10 percent (Drummond et al. 2014). For 4 Vesta, a reference spheroid with 
semiaxes ae = be ^ 285 km and Ce ~ 229 km provides an approximated description 
of the surface ( Janmann et al.[ |2012 ). These deviations from a spherical shape 
produce a difference of a few cm/s at most owing to the rather high inclination 
of the spin axes of 1 Ceres and 4 Vesta to the plane of the ecliptic (cf. Table [^. 
In the case of 2 Pallas, the approximating ellipsoid has semiaxes: a© = 275 ± 4, 
be = 258 ± 3, and Ce = 238 ± 3 km (Carry et al. 2010). Together with the low 


inclination of its spin axis to the plane of the ecliptic, this may produce a larger 
deviation from the results computed with a spherical shape, of the order of a few 
tens of cm/s. 

Other reflecting bodies often used to measure the solar RV are the Galilean 
satellites. Their parameters are listed in TableTheir spin axes are almost aligned 
with their orbital angular momenta - the maximum deviation is 0.46 degrees 
for Europa (e.g., Henrard fc Schwanen[ |2004[ ) - and their rotation periods are 
synchronized with their orbital periods owing to the strong tidal interaction with 
Jupiter. Since their orbits are in the equatorial plane of Jupiter that is inclined 
by less than ~ 4° 30^ to the ecliptic, their spin axes are almost orthogonal to the 
ecliptic plane and Z0 ~ zq ~ 90°. Moreover, their distance from the Sun is larger 
than in the case of the main-belt asteroids, thus Omax ~ 14°. As a consequence, 
their maximum radial velocity variation is found to be of a few cm/s. 


4 Discussion and conclusions 

We provide a method to compute the radial velocity perturbation induced by the 
rotation of a solar system body on the reflected solar spectrum. We focus on the 
case of a body with a uniform albedo and find that the perturbation depends in 
a complex way on the direction of its spin axis and the angle Sun-body-Earth. 
We treat both the case of a spherical reflecting body and of a triaxial ellipsoid. 
The introduced approach is based on vector scalar products rather than on the 
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Asteroid 

Omax 

(deg) 

ARVma,x 

(m/s) 

1 Ceres 

23.5 

1.31 X 10-4 

2 Pallas 

28.6 

2.428 

3 Juuo 

30.9 

0.361 

4 Vesta 

28.3 

0.060 

7 Iris 

28.9 

0.602 

20 Massalia 

29.6 

3.89 X 10-3 


Table 2 Maximum RV perturbation for the asteroids in Table II 


application of spherical trigonometry. With a suitable choice of the reference frame, 
we make the Cartesian components of the relevant vectors very simple so that the 
scalar products are easily computed. 

In the particular case when the angles z© and of the body spin axis to the 
directions of the Sun and the observer, respectively, are equal, we show that the 
radial velocity perturbation vanishes. On the other hand, when the inclination of 
the body spin axis to the plane of the ecliptic is less than about ±40° and the 
Sun-body-Earth angle a > 30°, the radial velocity perturbation can reach up to 



Fig. 5 Maximum values of the ratio ARV/Veq vs. the ecliptic latitude /3s of the spiu pole for 
a spherical asteroid. Differeut liuestyles iudicate differeut values of the augle a as labelled. 
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Satellite 

R 

(km) 

■Prot 

(day) 

Veq 

(m/s) 

lo 

1821 

1.77 

74.84 

Europa 

1605 

3.55 

31.98 

Ganimede 

2631 

7.16 

26.72 

Gallisto 

2411 

16.69 

10.50 


Table 3 Galileian satellite parameters. 


~ 0.1 of the equatorial rotation velocity of the body (cf. Fig. [^. In the case of 
asteroid 2 Pallas that has a spin inclination of ~ 30°, the perturbation can reach 
~ 2.4 m/s. For other bright asteroids, such as 1 Ceres and 4 Vesta, the effect is 
remarkably smaller, thanks to the rather high inclination of their spin axes to the 
plane of the ecliptic (cf. Table[^. However, in the case of 4 Vesta, it is necessary to 
take into account the present effect when an accuracy of a few cm/s is required. An 
easy estimate of the amplitude of the effect of rotation can be obtained from Fig.[^ 
when the inclination /3s of the spin axis to the ecliptic plane and the maximum 
value of the Sun-body-Earth angle Omax are known (cf. Sect.[^. 
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